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Abstract. We survey some recent developments in the theory of 
orthogonal polynomials defined by differential equations. The key 
finding is that there exist orthogonal polynomials defined by 2nd 
order differential equations that fall outside the classical families 
of Jacobi, Laguerre, and Hermite polynomials. Unlike the clas- 
sical families, these new examples, called exceptional orthogonal 
polynomials, feature non-standard polynomial flags; the lowest de- 
gree polynomial has degree m > 0. In this paper we review the 
classification of codimension m = 1 exceptional polynomials, and 
give a novel, compact proof of the fundamental classification the- 
orem for codimension 1 polynomial flags. As well, we describe the 
mechanism or rational factorizations of 2nd order operators as the 
analogue of the Darboux transformation in this context. We finish 
with the example of higher codimension generalization of Jacobi 
polynomials and perform the complete analysis of parameter values 
for which these families have non-singular weights. 



Even though the orthogonal polynomials of Hermite, Laguerre and 
Jacobi arose from various applications in applied mathematics and 
physics, these three families now serve as the foundational examples 
of orthogonal polynomial theory. As such, these classical examples 
admit many interesting and fruitful generalizations. 

A key property of classical orthogonal polynomials is the fact that 
they can be defined by means of a Sturm-Liouville problem. One of 
the foundational results in this Theorem of Bochner [5] which 

states that if an infinite sequence of polynomials of degree 0, 1,2, . . . 
satisfies a second order eigenvalue equation of the form 
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then p{x),q{x) and r{x) must be polynomials of degree at most 2,1 
and respectively. In addition, if the eigenpolynomial sequence is L^- 
orthogonal relative to a measure with finite moments, then it has to be, 
up to an affine transformation of x, one of the three classical classical 
families hsted above [B |29l [121 1271 [26] 

Implicit in the above definition of classical polynomials is the as- 
sumption that orthogonal polynomials form a basis for the standard 
polynomial flag C Pi C P2 C . . ., where Vn = span{l, x, x^, . . . , x"'}. 
In a pair of recent papers [lUllSIj, we showed that there exist complete 
orthogonal polynomial systems defined by Sturm-Liouville problems 
that extend beyond the above classical examples. What distinguishes 
our hypotheses from those made by Bochner is that the eigenfunction 
corresponding to the lowest eigenpolynomial of the sequence need not 
be of degree zero, even though the full set of eigenfunctions still forms 
a basis of the weighted space. 

Already some 20 years ago, particular examples of Hermite-like poly- 
nomials with non-standard flags were described in the context of su- 
persymmetric quantum mechanics [S1[S[. The last few years have seen 
a great deal of activity in the area of non-standard flags; the topic now 
goes under the rubric of exceptional orthogonal polynomials. There are 
applications to shape-invariant potentials [30], supersymmetric trans- 
formations [22], to discrete quantum mechanics [36], mass- dependent 
potentials [28], and to quasi-exact solvability [12]. As well, there are 
now examples of orthogonal polynomials that are solutions of 2nd order 
equations and that form flags of arbitrarily high codimension [35] . 

In light of the activity surrounding orthogonal polynomials with non- 
standard flags, we hope that it will be useful to summarize some key 
results and findings, and to supply stand-alone proofs to some key 
propositions. We note that the adjective "exceptional" was introduced 
in the context of our investigation of the equivalence problem for poly- 
nomial subspaces [201 H] • 



2. Preliminaries 

2.1. Polynomial flags. Let f/i C f/2 C . . ., where dimUk = k he a 
flag of real, finite-dimensional, polynomial subspaces. Let denote 
the degree of Uk] that is rik is the maximum of the degrees of the 
polynomials p E Uk- Let ik = nk + I — k denote the codimension of Uk 
in Vn^. , where the n + 1 dimensional polynomial vector spaces 

Vn = span{l, X, x^, . . . , x"} 
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make up the standard polynomial flag: Vq d Vi d 7^2 ■ ■ ■■ We say 
that {pk\'k=i is a basis of the flag if 

C/fc = span{pi,...,pfc} (2) 

Note that no generality is lost if we assume that degp^ — Uk- 

Definition 2.1. We call W = {Uk}kLi ^ degree-regular flag if < n^+i 
for all k. Equivalently, a flag is degree-regular if it admits a basis 
{PkjkLi such that degpfe < degp^+i for all k. 

Proposition 2.1. A polynomial flag is degree-regular if and only if the 
codimension sequence obeys £k < ik+i for all k. 

Henceforth, we will always assume that all flags under discussion are 
degree-regular. 

Definition 2.2. We say that the codimension of a polynomial flag is 
semi-stable if the codimension sequence ik admits an upper bound. 
If so, we call £ = limfc^fc the codimension of the flag. We say that 
codimension is stable ii ik — i is constant for all k. 

Likewise, whenever we speak of the codimension of a polynomial flag, 
we always assume that the flag has a semi-stable codimension. 

Let 1^2 (W) denote the vector space of 2nd order differential operators 
with rational real-valued coefficients 

T{y) = p{x)y" -\- q{x)y' -\- r{x)y, p,q,r&M.{x) (3) 

that preserve U; i.e., T{Uk) C Uk for all k. If p,q,r are polynomials, 
we call T a polynomial operator. Equivalently, an operator is non- 
polynomial if and only if it has a pole. 

Definition 2.3. We say that U is imprimitive if it admits a common 
factor; i.e., U is spanned by {qpk} where q{x) is a polynomial of degree 
> 1. Otherwise, we call U primitive. 

At this juncture, it is important to state the following two Propositions. 

Proposition 2.2. Let T{y) = py" -\- qy' -\- ry be a differential operator 
such that 

T{yi) = gi, i^ 1,2,3, 
where yi,gi are polynomials with yi, linearly independent. Then, 
p, q, r are rational functions with the Wronskian 

(y'i y[ yi\ 

W{yi,y2,y3) = det I y'^ y'2 2/2 

XVs y's yz) 

in the denominator. 
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Proof. It suffices to apply Cramer's rule to solve the linear system 



Vi Vi yi 



,11 



2/2 y2 1/2 

1/3 2/3 





□ 



We should also note that there is a natural linear isomorphism between 
the space 2^2 (W) of a flag U spanned by {pk}'k=i space V2{IA) 

of the imprimitive flag U spanned by {qpk}'k'=i- 

Proposition 2.3. LetT{y) be an operator that preserves U . Then, the 
gauge- equivalent operator T = qTq^^ preserves U. 

Proposition 12.21 makes clear why we restrict our definition of P2 
to operators with rational coefficients. Proposition 12.31 explains the 
need for primitive flags; these are the canonical representatives for the 
equivalence relation under gauge transformations. 

Definition 2.4. Let W be a polynomial flag of semi-stable codimension. 
We say that U is an exceptional flag if U is primitive and if 2^2 (W) does 
not preserve a polynomial flag of smaller codimension. 

Here are some examples to illustrate the above definitions. 

Example 2.1. The codimension 1 flag spanned by 1, a;^, x^, . . . , is ex- 
ceptional because the non-polynomial operator 

T{y) = y" - 2y'/x 

preserves the flag. Since £1 = 0,^^ = 1, k> 2, the codimension of this 
flag is not stable, but only semi-stable. 

Example 2.2. By contrast, the flag spanned by a; -|- l,a;^,x'^, . . . has 
a stable codimension 1 = 1. This flag is exceptional because it is 
preserved by the non-polynomial operator 

T{y) = y" -2{l + l/x)y' +{2/x)y. 

Example 2.3. Let hk{x) denote the degree k Hermite polynomial. The 
codimension 1 flag spanned by hi, h2, h^, . . . is not exceptional. The flag 
is preserved by the operator T{y) = y" — xy', but this operator also 
preserves the standard, codimension zero, polynomial flag. 

Example 2.4. The codimension 1 flag spanned by not 
exceptional because x is a common factor. The flag is preserved by the 
operator 

f{y) = y" - 2y'/x + 2y/x\ 
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However, observe that T = xTx ^ where T{y) = y" is an operator that 
preserves the standard, codimension zero, polynomial flag. 

Example 2.5. Let 

y2k-i = x''-' - i2k - y^k = - kx\ A; = 2, 3, 4, . . . (4) 

Consider the codimension-2 flag spanned by 1, 2/3, 1/4, 1/5, yg, ■ • ■• The 
degree sequence of the flag is 0,3,4,5, . . .; therefore the codimension 
is not stable, but merely semi-stable. The flag is preserved by the 
following operators |17j : 

TM = {^^ - 1)1/" - 2x1/', (5) 
T2(y)=x/-2(l + 2/(x2-l)y, (6) 
ri(y)=y" + x(l-4/(x2-l)y. (7) 

The flag is exceptional, because Ti and T2 do not preserve the standard, 
codimension zero flag. These operators cannot preserve a codimension 
1 flag, because, as will be shown in Lemma 13.21 an operator that pre- 
serves a codimension 1 flag can have at most 1 pole. 

2.2. Orthogonal polynomials. 

Definition 2.5. We will say that a 2nd order operator T{y) is exactly 
solvable by polynomials if the eigenvalue equation 

ny) = Xy. (8) 

has inflnitely many eigenpolynomial solutions y = yj with 

deg yj < deg yj+i, j = 1, 2, . . . . 

Let / = (xi,X2) be an open interval (bounded, unbounded, or semi- 
bounded) and let W{x)dx be a positive measure on I with flnite mo- 
ments of all orders. We say that a sequence of real polynomials {yj}j^i 
forms an orthogonal polynomial system (OPS for short) if the polyno- 
mials constitute an orthogonal basis of the Hilbert space L'^{I ,Wdx). 
If the codimension of the corresponding flag is stable, or semi-stable, 
we will say that the OPS has codimension m. 

The following deflnition encapsulates the notion of a system of or- 
thogonal polynomials deflned by a second-order differential equation. 
Consider the boundary value problem 

- {PyJ + Ry = XWy (9) 
lim (Py'u - Pu'y){x) = 0, z = 1, 2, (10) 
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where P{x), W{x) > on the interval / = (a;i, X2), and where u{x) is a 
fixed polynomial solution of We speak of a polynomial Sturm- 
Liouville problem (PSLP) if the resulting spectral problem is self- 
adjoint, pure-point and if all eigenfunctions are polynomial. 

If the eigenpolynomials span the standard fiag, then we recover the 
classical orthogonal polynomials, the totality of which is covered by 
Bochner's theorem. If the solution fiag has a codimension m > 0, 
Bochner's theorem no longer applies and we encounter a generalized 
class of polynomials; we name these exceptional, or polynomials. 

Given a PSLP, the operator 

T{y) = W-\Pyy-W-'Ry 

is exactly solvable by polynomials. Letting p{x),q{x),r{x) be the ra- 
tional coefficients of T{y) as in ([3]), we have 

P{x) = exp q/pdx^ , (11) 

W{x) = (P/pKx), (12) 
R{x) = -{rW){x), (13) 

Hence, for a PSLP, P{x), R{x), W{x) belong to the quasi-rational class 
[15], meaning that their logarithmic derivative is a rational function. 

Conversely, given an operator T{y) exactly solvable by polynomials 
and an interval / = {xi,X2) we formulate a PLSP (Q by employing 
( !TT|) -( [T3|) as definitions, and by adjoining the following assumptions: 

(PSLPl) P(x), W{x) are continuous and positive on I 
(PSLP2) Wdx has finite moments: Jj x^W{x)dx < 00, n>0 
(PSLP3) lim^^^^ P(x)x" = 0, i = l,2, n > 
(PSLP4) the eigenpolynomials of T{y) are dense in L^(/, Wdx). 

These definitions and assumptions (PSLPl), (PSLP2) imply Green's 
formula: 



T{f)g Wdx - r T{g)f Wdx = P{f'g - fg') (14) 



By (PSLP3) if f{x),g{x) are polynomials, then the right-hand side is 
zero. If / and g are eigenpolynomials of T{y) with unequal eigenval- 
ues, then necessarily, they are orthogonal in L'^{I,Wdx). Finally, by 
(PSLP4) the eigenpolynomials of T{y) are complete in L^(J, Wdx), and 
hence satisfy the definition of an OPS. 
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3. CODIMENSION 1 FLAGS 

The key result in this paper is the following. An analogous theorem 
for polynomial subspaces, rather than flags, was proved in [20] . 

Theorem 3.1. Up to affine transformations, the flag spanned by {x + 
l,x^,x^, . . .} is the unique stable codimension 1 exceptional flag. 

The proof proceeds by means of two lemmas. 

Lemma 3.1. ^4 primitive, codimension 1 polynomial flag is exceptional 
if and only if T>2 includes a non-polynomial operator. 

Proof. It's clear that a non-polynomial operator cannot preserve the 
standard polynomial flag. The standard flag is the unique codimension 
zero flag. Therefore, if a non-polynomial operator preserves a primitive, 
codimension 1 flag, that flag must be exceptional. 

Let us prove the converse. Let T[y) be a polynomial operator that 
preserves a polynomial flag. Consider the degree homogeneous decom- 
position 

N 

ny) = Yl T'^^y) (15) 

d=-2 

where 

Td{y) = x\adx^y" + M + 7d2/) (16) 

and where T/v is non-zero. Since the operator is polynomial, we must 
have (3-2 = 7-2 = 7-i = 0. Also note that 

= (J(J - + l^dj + 7d)a;^'+'- (17) 

Let 

yk{x) = x"'-' + lower deg. terms, k = 1,2, . . . 

be a basis of the flag. If > 0, then the leading term T^^y) raises 
degree, and hence 

TNix^') = 

for all k. By f ll7p . Ti^{y) annihilates at most 2 distinct monomials, 
a contradiction. Therefore, A^ < 0, and we conclude that T{y) is a 
Bochner-type operator ([1]). However, such an operator preserves the 
standard polynomial flag. Therefore, if U is exceptional, then there is 
at least one operator in V2(U) that doesn't preserve the standard flag. 
By the above argument, this operator must be non-polynomial. □ 

Lemma 3.2. Let T{y) be a non-polynomial operator that leaves invari- 
ant a codimension 1 polynomial subspace U cVn- If n > 5, then T{y) 
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has exactly one pole. Furthermore, up to a translation in x, a basis of 
U has one of the following three forms: 

1 + ax,x^,x^, . . . (19) 
l + aix^x + a2X^x^,...,x'' (20) 

Proof. By applying a translation if necessary, there is no loss of gener- 
ality in assuming that x = is a pole of T{y). Since the codimension 
is 1, the subspace admits an order- reduced basis of the form 

y. = + aj.ix", j = 1, . . . , z/ 

yj = x\ j = iy + l...n, 

where < u < n. The matrix representation of this basis is a n x n 
matrix in row-reduced echelon form. This matrix has n pivots and 1 
gap in position u. 

Our first claim is that u < 2. Suppose not. Then 

W{yi,y2, yz) = W{1, x, x^) + higher degree terms = 2 + 0{x). 

By assumption, T{yi) is a polynomial. Hence, Proposition 12.21 implies 
that a; = is not a pole of the operator, a contradiction. 

Having established that z/ < 2, we observe that X ^ Ob J Oil ^ G U. The 
Wronskian of these monomials is a multiple of x^ . Therefore, by Propo- 
sition 12.21 a; = is the unique pole. □ 

We are now ready to give the proof of Theorem 13.11 Let lA be an 
exceptional polynomial flag with stable codimension 1. By Lemma [STTj 
there exists a non-polynomial operator T{y) that preserves the flag. By 
Lemma [3^ this operator has a unique pole. Without loss of generality 
we assume that a; = is the unique pole. We rule out possibility f fTSj) . 
because if this holds for even one [7^, k > 5, then it must hold for all k. 
This would imply that X ^ Ob ^ OC ^ • • • IS 9j basis of the flag — a violation 
of the primitivity assumption. 

Let us rule out possibility ( 120|) . Let T{y) be a non-polynomial op- 
erator that preserves the flag. Since a: = is the unique pole, we can 
decompose the operator into degree-homogeneous terms 

ny) = X; Uy) 

d=Ni 

where Td{y) has the form shown in (fT6!) . and where T/v^ and Tjvj are 
non-zero. By the argument used in the proof of Lemma 13. H we must 
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have N2 < 0. By assumption, T(l + aix^),T{x + a2X^),T{x^) are 
polynomials. However, if Ni < —4, then this condition requires that 



which means that T/v^ is zero — a contradiction. If A^^i = —3, then T_3 
annihilates 1 contradiction. If A^^i = —2, then T_2 annihilates 

l,x,x^, another contradiction. Similarly, if A^^i = —1, then T_i annihi- 
lates 1. Hence T_i(|/) = a-ixy" + P-iy' . This means that there is no 
pole — a contradiction. Therefore, A^^i = A^2 = 0, but that means that 
T(|/) is a polynomial operator — again, a contradiction. This rules out 
possibility fljoj) . 

This leaves f|T9|) as the only possibility. Since we assume that the 
codimension is stable, deg(l + ax) = 1 and hence a 7^ 0. We scale x 
to transform 1 + ax to a; + 1. This concludes the proof of the main 
theorem. 

3.1. Xi polynomials. The above theorem explains the origin of the 
adjective "exceptional" and leads directly to a more general class of 
orthogonal polynomials — outside the class described by Bochner's 
theorem. Two families of orthogonal polynomials arise naturally when 
we consider codimension 1 flags. We describe these Xi polynomials 
below. 

In order to construct codimension 1 polynomial systems, we must 
consider V2 of the flag spanned by x + 1, x^, x^, . . .. 

Proposition 3.1. The most general 2nd order operator that preserves 
the flag {x + 1, x^, x^, . . .} has the form 



T{y) = {k2x''+k,x+k^)y"-{x+l) (k, + ^] y'+(k^ + ^] y, (21) 



where k2, fci, O'l"^ '^ea/ constants. 

See fl9[ Proposition 4.10] for the proof. 

Thus to obtain Xi orthogonal polynomials it suffices to determine 
all possible values of ^2,^1,^0 for which P{x),W{x) as given by ( ITTi) 
( IT^ satisfy the conditions of a PSLP. This analysis is performed in 
f2T\ . In summary, non-singular weights arise only for the case where 
/c2X^ -|- kix + ko either has two distinct real roots or if k2 = 0, fci 7^ 0. 
The first case leads to the Jacobi Xi polynomials; the second leads 
to the Xi Laguerre polynomials. In both cases, the polynomial fiags 
span a dense subspace of the respective Hilbert space [211 Proposition 
3.1, Proposition 3.3]. We summarize the key properties of these two 
families below. 



T^,(l) = T^,(x) =TjVi(x=^) =0, 
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3.2. Xi-Jacobi polynomials. Let a 7^ /3 be real parameters such that 
a > -1, /3 > -1, sgna = sgn/3. (22) 

Set 

a=i(/3-"), &=|^, c = h + l/a. (23) 
2 p — a 

Note that, with the above assumptions, we have |6| > 1. Now let 

Ml = X - c, Ui = {x- by, i > 2. (24) 
Define the measure dfia^i3 = Wa^p dx where 

W^c., = XG(-1,1). (25) 

[x — by 

Since W^,/? >Ofor— l<x<l, the scalar product 

if,g)a,p-= J fix)g{x)dfia,f}, (26) 

is positive definite. Also note that the above measure has finite mo- 
ments of all orders. 

We now define the Xi-Jacobi polynomials P^'^\ n = 1,2,... as 
the sequence obtained by orthogonalization of the fiag spanned by 
ui,U2,U3, . . . with respect to the scalar product (l26l) . and by impos- 
ing the normalization conditioi£| 

Pt'\l) = n["^l-y (27) 

From their definition it is obvious that degPi"'^'' = n. However, as 
opposed to the ordinary Jacobi polynomials, the sequence starts with 
a degree one polynomial. 

Next, define the 2nd order operator 

T^Ay) = - 1)1/" + 2a (\^) ((^ - c)y' - y) , (28) 

An elementary calculation shows that this operator preserves the fiag 
{ui} as defined above. Indeed, the fiag {ui} and the operator (1251) are 
obtained from the fiag x + l,x'^,x^, . . . and from the operator (121 p via 
the following specialization and affine transformation 

k2 = l, ki = -2ab, ko = {1 - b'^)a, x^-a{x-b). (29) 



^This convention differs from the one adopted in |21j. The change is made to con- 
form with the convention adopted below for the generahzed Xm Jacobi polynomials. 
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Multiplying both sides of the equation 

-TaAv) = 

by Wa^j3 leads to the following PSLP: 

((1 - x^)W^^p yj + 2a (^^j^^ W^,p 2/ = A W^^p y, (30) 
lim (l-x)"+^((x-c)y'-y) = 0, (31) 

a;— 

Urn {l + x f+\{x-c)y' -y) =0. (32) 

a;— 1+ 

The boundary conditions select polynomial solutions. The self-adjoint 
form of fl5U]) ensures that the solutions are orthogonal relative to dfla,i3- 
Therefore, the Xi Jacobi polynomials can also be described as polyno- 
mial solutions, y = P!f'^\x), of the following 2nd order equation: 

Ta,/3(y) = (n-l)(a + /3 + n)y. (33) 

3.3. Xi-Laguerre polynomials. For k > 0, set 

y^ = x + k + l, Vi = {x + k)\ i>2 (34) 
Define the measure dfik = Wk dx where 

Wu= ^^I'^ly . xG(0,oo) (35) 

and observe that Wk > on the domain in question. Therefore, the 
following inner product is positive definite: 

/>oo 

{f,9)k-= f{x)g{x)dfik, (36) 

^0 

Also note that the above measure has finite moments of all orders. 

We define the Xi-Laguerre polynomials L^^i, i = 1, 2, 3, ... as the se- 
quence obtained by orthogonalization of the flag spanned by f i, f2, fs, . . . 
with respect to the scalar product (!36l) and subject to the normalization 
condition 

Lk nix) = + lower order terms n > 1. (37) 

{n — 1)! 

Again, since we are orthogonalizing a non-standard flag, the Xi-Laguerre 
polynomial sequence starts with a polynomial of degree 1, rather than 
a polynomial of degree 0. 
Define the operator 

n{y) = -xy" + (l^) {{x + k + l)y' - y) (38) 
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and note that this operator leaves invariant the flag spanned by the f j. 
Indeed the flag {fj} and f l55]) are obtained from the flag x + . . . 
and (!2T]) via the specializations and an affine transformation shown 
below: 

^2 = 0, ki = —1, ko = k, X I— )■ X + /c (39) 
The corresponding PSLP takes the form 

{xWkyJ + (l^) WkV = Xy, (40) 
\im x''^\{x + k+l)y' -y) =0, (41) 
lim e-^((x + A; + l)?/'-?/) = 0. (42) 

a;— ^oo 

As before, the boundary conditions select polynomial solutions, while 
the self-adjoint form of ( HOjl ensures that these solutions are orthogonal 
relative to dfik- Therefore, the Xi Laguerre polynomials can also be 
defined as polynomial solutions, y = Lk^n, of the following 2nd order 
equation: 

Tk{y) = {n- l)y. (43) 

Having introduced the Xi Jacobi and Laguerre polynomials, we are 
able to state the following corollary of Theorem 13.11 The proof is 
found in \21\ . 

Theorem 3.2. The Xi- Jacobi polynomials and the Xi Laguerre poly- 
nomials are the unique orthogonal polynomial systems defined by a sta- 
ble codimension- 1 PSLP. 



4. Higher codimension flags 

Even though the first examples of exceptional orthogonal polynomi- 
als involve codimension- 1 flags, recently announced examples [35] are 
proof that exceptional orthogonal polynomials can span flags of arbi- 
trarily high codimension. Another important development is the recent 
proof [30] that the codimension- 1 families can be related to the clas- 
sical orthogonal polynomials by means of a Darboux transformation. 
In a follow-up publication [22j, it was shown that the new higher codi- 
mension examples are systematically derivable by means of algebraic 
Darboux transformations [17j . 

A closely related development involves the notion of shape-invariance 
[T3], a methodology related to the study of exactly solvable poten- 
tials. The close connection between solvable potentials and orthogonal 
polynomials is well recognized; consider the relationship between the 
harmonic oscillator and Hermite polynomials, for example. 
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In the orthogonal polynomial context, we have to factorize general 
second order operators, not just Schrodinger operators. The shape- 
invariance property of the classical potentials manifests as the Ro- 
drigues' formula for the corresponding polynomials. It has been pointed 
out that all the potentials related to exceptional orthogonal polynomi- 
als exhibit the shape-invariance property |30[ [38] , and that therefore, 
like their classical counterparts, exceptional orthogonal polynomials 
have a Rodrigues' formula. This phenomenon was studied and ex- 
plained in ^J, where it was shown that the shape-invariance property 
of the Xjji (exceptional, codimension m) polynomials follows from the 
permutability property of higher-order Darboux transformations. 

4.1. The Darboux transformation. In the remainder of this sec- 
tion, we review the some key definitions and results related to Darboux 
transformations and shape-invariance, and then illustrate these ideas 
with the example of Xm Jacobi polynomials [39]. 
Consider the differential operators: 

T{y) = py" + qy' + ry (44) 

A{y) = b{y' -wy), B{y) = b{y' - wy), (45) 

where p, q, r,b,w,b,w are rational functions. 

Definition 4.1. We speak of a rational factorization if there exists a 
constant Aq such that 

T = BA + Xo (46) 
If the above equation holds, we call 

f = AB + Aq. (47) 
the partner operator. We call 

0(x) = exp J wdx, w = (j)'/(j) (48) 

a quasi-rational factorization eigenfunction and b{x) the factorization 
gauge. 

The reason for the above terminology is as follows. By (H6l) . 

T(0) = Ao0; (49) 

hence the term factorization eigenfunction. Next, consider two fac- 
torization gauges bi{x),b2{x) and let Ti{y),T2{y) be the corresponding 
partner operators. Then, 

T2 = /i~^Ti/x, where fi{x) = bi{x) /b2{x). 
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Therefore, the choice of b{x) determines the gauge of the partner op- 
erator. This is why we refer to b{x) as the factorization gauge. Also, 
note that in [17] the above construction was referred to as an algebraic 
Darboux transformation. However, in hght of the recently recognized 
role played by operators with rational coefficients, the term rational 
factorization seems to be preferable. 

Proposition 4.1. Let T{y) be a 2nd order operator exactly solvable by 
polynomials, and let(j){x) be a quasi-rational factorization eigenfunction 
with eigenvalue Aq • Then, there exists a rational factorization fl46p such 
that the partner operator is also exactly solvable by polynomials, and 
such that the partner flag is primitive (no common factors) . 

Proof. Let w{x) = cj)' (x) / (f){x) and let b{x) be an as yet unspecified 
rational function. Set 

w = —w — q/p + b' /b, (50) 
b = p/b, (51) 

and let A{y),B{y) be as shown in in (H5|l . An elementary calculation 
shows that (l46l) holds. Let yi,y2, ■ ■ ■ be a degree-regular basis of the 
eigenpolynomials of T{y). We require that the flag spanned by A{yj) 
be polynomial and primitive (no common factors). Observe that if 
we take b{x) to be the reduced denominator of w{x), then A{yj) is a 
polynomial for all j. However, this does not guarantee that A{yj) is 
free of a common factor. That is a stronger condition, one that flxes 
b{x) up to a choice of scalar multiple. Finally, the intertwining relation 

TA = AT (52) 

implies that the A{yj) are eigenpolynomials of the partner T. □ 

Finally, let us derive the formula for the partner weight function. 

Proposition 4.2. Suppose that a PSLP operators T{y) = py"+qy'+ry 
is related to a PSLP operator T{y) = py" + qy' + fy by a rational fac- 
torization with factorization eigenfunction (f){x) and factorization gauge 
b{x). Then the dual factorization gauge, factorization eigenfunction and 
weight function are given by 

bb = p (53) 
W/b = W/b, (54) 
60 = l/{W(j)) (55) 
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Proof. Equation (153|) follows immediately from ( H5|) ( l46l) . Writing 

^'(i/) = py" + Qy' + ry, (56) 
equation ( HTll implies that 

w + i& = — g/p + = —q/p + &'/^- (5''') 

Hence, 

q = q + p' - 2pb'/b. (58) 

From here, ( l54l) follows by equations ( ITTl) ( |T2l) . Equation ( l55i) follows 

from (Ml. □ 



The adjoint relation between A and -B allows us to compare the 
norms of the two families. 

Proposition 4.3. Let T{y),T{y) be PSLP operators related by a ra- 
tional factorization fj46|) P7j) . Let {yj} be the eigenpolynomials ofT{y) 
and letyj = A{yj) be the corresponding partner eigenpolynomials. Then 

2 TJ/ J™ I \ \ \ „,2 



xi J x\ 



A{yf)' Wdx = (Ao - A^) / y; Wdx (59) 



where Xi < x^ are the end points of the Sturm- Liouville problem in 
question. 

Proof. As a consequence of fl54p . A and —B are formally adjoint relative 
to the respective measures: 



X2 rx2 



X2 



A{f)g Wdx + / B{g)f Wdx = {P/b)fg , (60) 

Xl J X\ ^'1 

where P{x) is defined by (fTTl) and where b{x) is the factorization gauge. 
By assumption both W{x),W{x) are positive on (xi,a;2). By flM|) . 
W = P/b"^. In particular, the numerator b{x) cannot have any zeroes 
in (xi, 0:2). As well, if either xi, X2 are finite, we cannot have b{xi) = 0, 
because that would imply that W is not square integrable near x = Xi. 
Hence |l/6(x)| is bounded from above on {xi,X2). Therefore, if f,g 
are polynomials then the right-hand side of flBUj) vanishes by (PSLP3). 
Therefore, 

fX2 rX2 fX2 



i>X2 rX2 rX2 

/ A{yjf Wdx = - B{A{y,))yj Wdx = (Aq - A,) / y] Wdx 

J Xl J Xl J XI 



(61) 

□ 
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4.2. Shape-invariance. Parallel to the spectral theory [21110], ra- 
tional factorizations of a solvable operator T{y) can be categorized as 
formally state-deleting, formally state-adding, or formally isospectral. 
The connection between these formal, algebraic Darboux transforma- 
tions and their analogues is discussed in [22j . We speak of a formally 
state-deleting transformation if the factorization eigenfunction is 
the lowest degree eigenpolynomial of T{y). We speak of a formally 
state-adding transformation if the partner factorization eigenfunction 



(x) = exp w dx 



with w defined by (l50ll . is a polynomial. We speak of a formally isospec- 
tral transformation if neither nor are polynomials. Examples of all 
three types of factorizations will be given below. 

State-adding and state-deleting factorizations are dual notions, in 
the sense that if the factorization of T is state-deleting, then the fac- 
torization of T is state-adding, and vice versa. 

Definition 4.2. Let n E K he a, parameter index set and let 

T^iy) = Pix)y" + q^{x)y' + r^{x)y, k e K, (62) 

be a family of operators that are exactly solvable by polynomials. If 
this family is closed with respect to the state-deleting transformation, 
we speak of shape-invariant operators. 

To be more precise, let tc^^x) = yK,i{x) be be the corresponding 
eigenpolynomial of lowest degree. Without loss of generality, we assume 
that the corresponding spectral value is zero, and let 

= B^A,, A^n^ = (63) 

be the corresponding factorization. Shape-invariance means that there 
exists a one-to-one map h : K K and real constants such that 

T;,(«) = A^B^ + A,. (64) 

In accordance with (fTTj) . define 

Pt,{x) = exp g^/p^ (65) 

Let &k(x) denote the shape-invariant factorization gauge; i.e.; 

A(2/) = (&KA«)>V(7r,,2/), (66) 

where 

mf,9) = f9'-f9 (67) 
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p(«,/3)ni 



denotes the Wronskian operator. By equation ( !64l) . 

qh(K) = Qk+p' - '2pb'^/b^. (68) 

It follows that 

pPjPhi.) = bl. (69) 

Below, we will use this necessary condition to derive the factorization 
gauge of a shape- invariant factorization. 

4.3. Jacobi polynomials. The operator 

To^M ■= (1 - ^')y" +(/?-« + (« + /? + 2)x)y', (70) 
preserves the standard flag, and hence is exactly solvable by polyno- 
mials. The classical Jacobi polynomials P^'^\x), a,/? > — 1, n = 
0,1,2,... are the corresponding eigenpolynomials: 

Ta.,pPlC'^'^ = -n{n + « + /? + l)Pt^\ (71) 
subject to the normalization 

+ a" 
n 

The orthogonality is relative to the measure Wa,(5{x)dx where 

WaA^) = (1 - a;)"(l + xf, X G (-1, 1). (72) 

The classical operators are shape-invariant by virtue of the following 
factorizations: 

Ta,l3 = Ba^isAa^is, (73) 

= Aa,f3Ba^i3 + a + (3 + 2 where (74) 
B^^y) = (1 - x^)y' + (/?-« + (a + + 2)x)y, (75) 
= (1 - x)-°(l + x)-^ {y{l - x)"+^(l + xf+^)' (76) 

A^Ay) = y'- (77) 

As a consequence, Ba^y) acts as a raising operator: 

B^^pPlC^^'^^'^ = -2(n + l)Pif,\ (78) 
and Aa^js^y) as a lowering operator: 

Pi"'^'^ ' = ^(1 + « + /? + n)P,+i,^+i,„_i. (79) 
The classical Rodrigues' formula, namely 

(1 _ _ ((1 _ _ 3;)/^+") = (-2)"n!p("'^)(x), 

(80) 
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follows by applying n iterations of the raising operator to the constant 
function. 

The quasi-rational solutions of Ta^p{y) = Xy are known Section 
2.2]: 

Mx) = P!n'^\x), Xo = -m{l + a + P + m) 

(j)2{x) = (1 - xy"{l + x)"^P_c.,-/3,m(a;) Ao = (1 + m)(a + /3 - m) 

03(x) = (1 - Xy"'P^a,l3,mix) Aq = (1 + /? + m) (« - m) 

04(x) = (1 - xy^Pa-f}^m{x), Aq = (1 + « + m) (/3 - m) 

where m = 0, 1, 2, . . .. The corresponding factorizations were analyzed 
in [16]. Of these, 0i with m = corresponds to a state-deleting trans- 
formation and underlies the shape-invariance of the classical Laguerre 
operator and the corresponding Rodrigues formula. For m > 0, the 
01 factorization eigenfunctions yield singular operators and hence do 
not yield novel orthogonal polynomials. The 02 family results in a 
state-adding transformation. The resulting flags are semi-stable, like 
in Example 12. 5j see [H] for a discussion. The type 03 04 factorizations 
result in novel orthogonal polynomials, provided a, /3 satisfy certain 
inequalities. These families were referred to as the Jl, J2 Jacobi poly- 
nomials in [37j. The two families are related by the transformations 
a ■v^ P,x ^ —X. We therefore focus only on the 03 factorization; no 
generality is lost. 

The derivations that follow depend in an elementary fashion on the 
following well-known identities of the Jacobi polynomials. We will ap- 
ply them below without further comment. 

P^(-'«(x) = l, p(-'P\x) = 0, n<-l, (81) 
P("'«(x) = (-irp(°'^)(-x), (82) 
{x - l)p(r'« '(x) = (a + m)Pj:~'''+'\x) - aPt'^\x) (83) 

pIC'"^ '(x) = i(l + a + /3 + n)P,+i,^+i,„_i(x), (84) 

Pt^'-'\x) - P^--'''\x) = Ptf\x) (85) 

4.4. The Xm Jacobi polynomials. Fix an integer m > 1 and a, /3 > 

— 1, and set 



Up,m = PL'''^\x) (86) 
Take 03 ( the factorization eigenfunction and take 

h{x) = (1 -X)^a,f3,m 
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as the factorization gauge. Applying ( H5i) ( l50i) and the identities ( IHTl) - 
dHS]) , we obtain the following rational factorization of the Jacobi op- 
erator ([28]): 

Ta,i3 = Ba,i3,mAa,i3,m - {m - a){m + (3 + 1) where (87) 

Aa,p,m{y) = (1 - X)^a,fi,m v' + {m - a)^a+l,fi+l,m V (88) 
Ba,^,m{y) = ((1 + ^)?/' + (1 + (§9) 

By fl5^ . 0(x) = (1 + a;)^^^'^ is the dual factorization eigenfunction. 
Since neither 03 (x) nor 0(a;) is a polynomial, ( 187|) is an example of a 
formally isospectral factorization. The corresponding partner operator 
is shown below: 

Ta,l3,m = Aa+l,(3-l,mBa+l,t3-l,m - (m - a - 1) (m + /3) , (90) 

Ta,i5,m{y) = TaAv) " 2pa+i,/3-i,m((l - x^)y' + 6(1 - x) y) (91) 
+ m(a — h — m + l)y, where 

Pa,l3,m = ia,l3,mlia,P,rn (92) 
= ^(l-« + /3 + H ^Q_l,/3+l,m--l/^Q,/3,m (93) 

Based on the above factorization, we define the X^, exceptional 
Jacobi polynomials to be 

( ^ \m+l 

piaM ^ ^_^L^ ^ pia+i,p-i) ^ J = ri-m>0 (94) 

a + 1 + J ■' 

1 + a + /3 + J p(a+2,/3) 



X - l)Pi-"-i''^-i)i^(!f (95) 



2(a + l+j) 

(96) 



I p{-2-a,l3)p{ 

a+l+j " ' 
By construction, these polynomials satisfy 

T„,/3,r„Pi"'^'™^ = -(n-m)(l + a + /3 + n- m)P^"''''"^) (97) 

With the above definition, the generalized Jacobi polynomials obey the 
following normalization condition 

p(-/'-)(l)=f" + "-"^V^Y n>m. (98) 



Note that the operators and polynomials extend the classic family: 

TaMy) = TaAv) (99) 
pia,m = p^",/3). (100) 
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The norms of the classical polynomials are given by 



-1 



where 



2"+/^+ir(a + 1 + n)r(/3 + 1 + n) 



Wo.,^A^) = ^; ^ ^7"^ (101) 



n\{a + l3 + 2n + l)T{a + /3 + n + 1) 
By (IMl) . the weight for the Jacobi polynomials is given by 

[1 -x)''{l + x)^ 

In order for orthogonality to hold for the generalized polynomials, 
we restrict a, (3 so that the denominator in the above weight is non-zero 
for — 1 < X < 1. We also want to avoid the degenerate cases where 
X = ±1 is a root of ^a+i,i3-i,m- To ensure that we obtain a codimension 
m flag, we also demand that deg,^Q,+i,^_i,m = The proof of the 
following Proposition follows from the analysis in jUl Chapter 6.72]. 

Proposition 4.4. Suppose that a,(3 > —1. Then deg^Q,+i,^_i,m = m 
and ^Q,+i,/3_i,m(±l) 7^ z/ and only if f3 and 

a,a-/3-m + 1 ^ {0,l,...,m- 1} (102) 

Proposition 4.5. Suppose that the conditions of the preceding Propo- 
sition hold. The polynomial ^a+i,p-i,m{x) has no zeros in (—1, 1) if and 
only if a > m — 2 and 

sgn(a — m + 1) = sgn(/3). 

This is a good place to compare the above results to the parameter 
inequalities imposed in [371 EH] ■ These references impose the condition 



a > (3 > m - 1/2. 



Unlike Proposition 14.51 this condition fails to describe the most general 
non-singular weight Wa^p^m- Consider the following examples: 

1^1/3,-1/2,2 - 288 (7^2 + 2x-41)^ 

W^u 1/2 2 = 128^ \ , ' — ^ \, 104 
■ (5x2 - 14x + 29)2 V ; 

Neither of the above examples satisfy the parameter inequalities of 
139] , but both weights are non-singular on (—1, 1) and have finite 
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moments of all orders. On the other hand, the parameter values m = 
2,a = 3/2, /3 = 1/2 give 

W^3/2,l/2,2(x) = ^(l-x)3/2(l + x)l/2, 
p(3/2,l/2,2) ^ 3^(3/2,1/2) ^ ^ q 

In other words, for certain singular values of the parameters, the codi- 
mension is actually less than m, and in some instances (such as the one 
above) even yield the classical polynomials. The condition f ll02p must 
be imposed in order to avoid such singular possibilities. 

Proposition 4.6. Suppose that a,(3 > —1 satisfy the conditions of 
Propositions 4-4 and \4-5\ The norms of the Xm Jacobi polynomials 
are given by 



^'u/ , {l + a + k-m){/3 + m + k) 



p{a,/3,m)/ X 



Proof. This follows directly from ( l59l) . □ 

We summarize the above findings as follows. 

Theorem 4.1. Let m > 1 and a, (3 > —1 be such that a > m — 
2, sgn(a — m + 1) = sgn(/3) and such that fll02p holds. Let U be 
the stable, codimension m flag spanned by polynomials y{x) such that 

{l + x)y' + /3y is divisible by Pm ^''^ ^\^)- Let Wa,i3,m{x) be the weight 
defined by fl5^ . Then, the Jacobi polynomials, as defined by flMl) 
are the orthogonal polynomials obtained by orthogonalizing the flag U 
relative to the weight Wa^fi^mix) and subject to the normalization con- 
dition 



Finally, let us discuss shape-invariance of the generalized Jacobi op- 
erators. The following Proposition was proved in [22] . 

Proposition 4.7. Let A,B be the operators defined in fllOSp fll09p . 

Then 

Ta,j3,m = -Ba,/3,m^Q,/3,?TU Aa,l3,mPm = (105) 

is the state- deleting factorization of the X^ Jacobi operator. Further- 
more, 

+ a + (3 + 2, (106) 

is the dual state-adding factorization. 

In essence, we are asserting that the generalized operators obey the 
same shape-invariance relations as their classical counterparts; c.f., 
equations (JTS]) ([74]). 
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The proof of Proposition 14.71 rehes on the permutabihty property 
of higher order Darboux transformation and goes beyond the scope of 
this survey. We hmit ourselves to exphcitly deriving the raising and 
lowering operators used in the above factorization. 

We already know the factorization eigenf unction: 

We make use of (!69l) to determine the factorization gauge. Making use 
of the fact that h{a, /3) = (a + 1, /3 + 1) we obtain 

K^) = 7 • (107) 

Sa+l,;S-l,m 

We use Proposition 14.21 to derive the dual factorization gauge 

b{x) = [l-X )— 

and the dual factorization eigenfunction, 

0(x) = (1 - a;)-'^-^(l + x)-^-iea+i,/3-i,m. 
We thereby obtain 

Aa,(3,miy) = /""^^'^'"^ iv' - Pa+2,l3,my) (108) 
'ia+l,l3-l,m 



Ba,l3,m{y) = {I- X 



,m 



Q!+2,/3,m 



, , a + l /3 + l\ 

(109) 



These shape- invariant factorizations serve as a good illustration of 
the duality between formal state-adding and state-deleting transfor- 
mations; here 0(x) is a polynomial but 0(x) is merely a quasi-rational 
function. 

As well, the shape- invariant factorization illustrates that 6(x), the 
factorization gauge, is not necessarily a polynomial. Here, 

W = -(j)'/(f) = Pa+2,f),m- 

The denominator is C,a+2,i3,m but the transformation 

y ^ ^a+2,p,Uy' - wy)^ y = ^i"'''''"^ ^ > ^ 

would produce an imprimitive flag. The common factor is ^a+i,ii-i,m, 
and that is why the correct factorization gauge is the rational function 
shown in (11071) . 
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As a consequence of these shape-invariant factorizations, we have the 
following lowering and raising identities for the Jacobi operators; 

c.f., m dzi 
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